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Abstract. This paper investigates questions related to the modularity 
in discrete models of biological interaction networks. We develop a the- 
oretical framework based on the analysis of their asymptotic dynamics. 
More precisely, we exhibit formal conditions under which agents of inter- 
action networks can be grouped into modules. As a main result, we show 
that the usual decomposition in strongly connected components fulfils 
the conditions of being a modular organisation. Furthermore, we point 
out that our framework enables a finer analysis providing a decomposi- 
tion in elementary modules. 

Keywords: modularity, interaction networks, discrete dynamics, equilib- 
ria. 



1 Introduction 

Understanding and exhibiting the relations between phenotypes and in- 
teractions in biological networks, i.e., the links between structures and 

functions [1], are among the most challenging problems at the frontier of 
theoretical computer science and biology. Many phenotypes can be asso- 
ciated to interactions of biological agents {e.g., genes or proteins) working 
together to guarantee some specific functions. This leads to group agents 
into modules and to associate to them one or more biological functions. It 
follows that biological networks can be seen, at a more abstract level, as 
modular networks in which interacting elements are modules that carry 
the biological information necessary to translate into precise functions. 

Modularity is present in various kinds of networks including metabolic or 
signalling pathways and genetic or protein interaction networks. Modular 
organisations are notably emphasised in embryogenesis [2,3] where mod- 
ules of genes are coordinated in the development process. Furthermore, 
methods related to modules discovery in interaction networks are gener- 
ally based on both the analysis of the networks structures (a field close 



* Corresponding author: franck.delaplace@ibisc.univ-evry.fr. 



to graph theory) and the study of their associated dynamics [4]. Struc- 
tural analysis identifies sub-networks with specific topological properties 
motivated either by a correspondence between topology and functional- 
ity [5,6] or by the existence of statistical biases with respect to random 
networks [7] . Specific topologies like cliques [8] , or more generally strongly 
connected components (SCCs) are commonly used to reveal modules by 
structural analysis. Particular motifs [9,10] may also be interpreted as 
modules viewed as basic components. They represent over-expressed bi- 
ological sub-networks with respect to random ones. Dynamical analysis 
lays on the hypothesis that expression profiles provide insights on the 
relationships between regulators, modules being possibly revealed from 
correlations between the biological agents expressions. For instance, us- 
ing yeast gene expression data, the authors of [11,12] inferred models of 
co-regulated genes and the condition under which the regulation occurs. 
As a consequence, the discovery of a modular organisation in biological 
interaction networks is closely related to the influence of the expression 
of agents on the others and needs to investigate their expression dynam- 
ics [13,14,15,16,17]. 

Using a discrete model of biological interaction networks [18,19], we pro- 
pose an approach that analyses the conditions of modules formation and 
characterises the relations between the global behaviour of a network and 
the local behaviours of its components. The chosen model of interaction 
networks is based on the assumption that phenotypcs arc related, at the 
molecular level, to characteristic states of some biological agents assimi- 
lated to equilibria^. Thus, using dynamical analysis, we show conditions 
under which interaction networks can be divided into modules so that the 
composition of these modules behaviours matches the global behaviours 
of the networks. 

The paper is structured as follows: First, Section 2 introduces the main 
definitions and notations used throughout the paper. Then, Section 3 
presents the central notion of a modular organisation of a network along 
with its structural and dynamical properties. Section 4 defines elemen- 
tary modular organisation and the conditions leading to obtain it. Some 
concluding remarks and perspectives are provided in Section 5. 



^ In systems biology, equilibria are usually associated to differentiated cellular 
types [20] and specific behavioural properties like lysis and lysogeny in the bac- 
teriophage A [21]. 
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2 Preliminaries 



First, we introduce basic notations. Let C S x S he & binary relation 
on S. Given s,s' G S and S' C S, we denote by s ^ s' the fact that 
(s, s') G by (s = {s' I s s'} the image of .s by and by 

(S" ^) its generalisation to S' . Similarly, we denote by s) and {-^ S') 
the corresponding preimages. The composition of two binary relations will 
be denoted by ^ o the reflexive and transitive closure by -^*. 

2.1 Interaction network and its associated dynamics 

Let us assume that a network 77 is composed of a set ^ = {ai, . . . ,a„} 
of agents. Each ai £ A has a local state, denoted by Sa., taking values in 
some nonempty finite set Sa^- A state (or a configuration) of ry is defined 
as a vector s G 5 associating to each £ A a value in 5a., where S = 
Sai X ... X Sa„ is the set of all possible states of the network. For any 
X C A and s G S, wc denote by s\-^ the restriction of s to the agents in 
X; this notation naturally extends to sets of states. 

An evolution of r/ is a relation C. S x S where each s ^ s' is a transition 
meaning that s evolves to s' by Thus, the global evolution of rj can be 
represented by a directed graph Q = (5, ^) called the state graph (or the 
transition graph). In this work, we pay particular attention to the notion 
of local evolution, since each agent a £ A has its own evolution The 
collection of all these local evolutions results in the asynchronous view of 
the global evolution of rj. 

Definition 1. The asynchronous dynamics (or dynamics for short) of 
a network rj is defined as the triple {A,S, (^a)aeA); where A is a set of 
agents, S is a set of states, and for each a £ A, S x S is a total or 

empty relation characterising the evolution of agent a such that for any 
s s' , either s = s' or s differs from s' only on the a-th component. 

We are now in a position to introduce formally the interaction network as 
a family of functions^ ij = {r]a}aeA, such that each r]a ■ S ^ Sa defines the 
next state r)a{s) with respect to the asynchronous evolution of a from s. 
Network rj allows to deduce a directed graph of interactions G = {A, — >) 
such that — > aj if occurs in the definition of ijay Notice that 77 
implies G, but the opposite is in general not true. When is empty 
for some a {i.e., the local state of a remains invariant), then a plays the 

* Actually, v ■ A ^ S ^ UaeA ^a- 
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role of an input, which means that no other agents of A influences it (i.e., 
there are no arcs towards a in G). 

Given a set S' C S, we introduce the notions of orbits, equihbria and 
attractors: 

— An orbit of S' is the set of states f2{S') = {S' -^*) comprising S' 
and all the states reachable from S"; 

— An equilibrium e G S* is a state endlessly reachable by ^{S') denotes 
the set of equilibria reachable from S'; 

— An attractor is a set of equilibria ECS such that e E: ^{{e}) = 
E. In a state graph, attractors are sets of states belonging to terminal 
strongly connected components that can be of two kinds: 

- a stable state is a singleton E C. S; 

- a limit set is an attractor E such that \E\ > 1. 

The restrictions of ^ to X C ^4, denoted by ^x= UaeX * 
state s, denoted by s\x, lead to the following definitions. 

Definition 2. Two states si,S2 G S are said to be X-equivalent and are 
denoted by si ~x ^2; for some X C A, if si\-^ = S2\x> if they are 
the same for all agents in X. 

This equivalence relation naturally extends to sets of agents, i.e., Si ~x 



Definition 3. We consider the following operators with the same signa- 
ture V{A) {ViS) V{S)), where X C A and S' C S: 

— the orbit operator f2, defined as Qx{S') = (5" —^*x), 

— the equilibria operator "F, defined as <Px (S") = {s G f2x{S') \ Vs' G 



Proposition 1 below emphasises specific properties of operator ^. 

Proposition 1. Let X C A be a subset of agents, has the following 
properties, for all sets S',S" subsets of states: 

a. Idempotency: ^x ° 'I'xiS') = ^x{^x{S')) = '^x(5"); 

b. Upper- continuity: ^x{S' U S") = ^x{S') U "^xiS"); 

c. Monotony (order-preserving) : S' C S" =^ <Px(S'') C ^x{S"). 

Proof. Let Eqx{s) be the predicate meaning that s is an equilibrium for 
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Fig. 1. An interaction network and its state graph. 



a. By expanding ^x{^x{S')), we have: 

Wx{^x{S')) = {se f}x{^x{S')) I Eqx{s)} 

= ^x{S'). 

b. By definition, ^x{S' U S") = {s e nx{S' U S") \ Eqx{s)}, where 
f2x{S' U S") = (S" U S") Since is upper-continuous on the 
lattice of state sets, we have: 

^x{S' U S") = {se f2x{S') U nxiS") I Eqxis)} 

= {s e nxis') 1 u {s G nxis") I 

= lFx(-S')U!f'x(5")- 

c. An upper-continuous function is monotonous. □ 

Example in Figure 1 illustrates the dynamics of a Boolean network^, i.e., 
in which for all a e A, Sa = {0, 1}. It is defined by rj, each r]a being the 
local transition function of agent a. Given a state s, the evolution s — 
means that s' is obtained by applying rja to s. More formally: 

s^as' = s'^ = riais) A [W eA\{a}: s'^ = s„) . 

Figure 1 (top) provides 77 and its graphical representation. Figure 1 (bot- 
tom) depicts the state graph Q of rj. More precisely, Q represents the 
dynamics of ry for each state s G {0, 1}^. Let us remark that: 

^ In the state graphs of interaction networks cited as examples, self loops are omitted. 
Stable states are depicted in grey while limit sets are in black. 
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Fig. 2. An interaction network, its state graph and the corresponding regulation graph. 

- the orbit of {1111} is I?({1111}) = {1111, 1101, 1100}; 

- 1100 is an equilibrium, as well as 0101 and 0010 are; 

- the set of equilibria reachable from 0000 is {Oxyz \ x,y,z e {0, 1}}; 

- two attractors exist: a stable state {1100} and a limit set tZ'({0000}). 

2.2 Regulation 

The regulation relation specifies a state-based dependence between two 
agents. Agent is a regulator of agent a^, denoted by — > a£, if at 
least one modification of the state of requires a modification of the 
state of ajfc. 

Definition 4. ak regulates ai, denoted by — > ae, if and only if there 

exist two states s, s' G S such that (s ~A\{afe} ^ag) (s' ^a^))- 

By extension, given Xi,Xj C A, Xi — > Xj if and only ifBa^ € Xi,3a^ G 
Xj : Ok — > a^. 

In Figure 1 , the sets of regulators of agents ai , 02, as and 04 are respec- 
tively {ai}, {01,03}, {02} and {03}. Notice also that there are the fol- 
lowing relations on sets of agents: {ai, 02} — > {as, 04}, {03} — > {02, 04} 
and {fli} — ?• {02,03} — 7- {04}. Another example of a regulation graph 
is given in Figure 2 (right). It shows that interaction (01,02) in the inter- 
action network is actually not a regulation because no modification of Oi 
influences the state of 02- All the other interactions are effective, meaning 
that the underlying regulation graph contains all interactions from the 
network but (01,02). 

3 Composition of equilibria 

In this section, we give rise to the relation between the equilibria of an 
interaction network and that of its parts. This allows to consider a modu- 
lar view of the system in which each part is assimilated to a module, i.e., 
a subset of agents with their dynamics. It means that modules, which in- 
fluence each other, reveal the underlying biological functions materialised 
by their equilibria. 
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IVai (s) = -'Sai V Sa2 
Va3 (s) = Sa2 V -iSag 




Fig. 3. An interaction network and its state graph. 



3.1 Modular organisation 

Our objective is to find a decomposition of A into a set of modules, i.e., a 
partition^ of A, together with a composition operator of the equihbria 
of these modules, allowing to retrieve the global equilibria of the network. 
It is actually a challenge to find an adequate operator 0, i.e., such that a 
decomposition into a set of modules {Xi, . . . , X^} satisfies: 

^x,(2)...(Z>^Xr^=^U^^^Xr (1) 

One can easily see that, in general, taking = U is not a solution. For 
example, if we consider the interaction network in Figure 3 and a partition 
into two modules {01,02} and {as}, then the corresponding sets of equi- 
libria are respectively {110,111} and {000,001,100,101,011,111}, while 
the set of global equihbria is {111}, which is not the union of the previous 
ones. However, one may see that, for the same modules, the computation 
of the equilibria of {0.3} from the equilibria of {04,02}, gives the desired 
pmpeity ^{a3}°'^{aua2} = '^{ai.aj.ag}' whereaS 'i'{ar,a2}°'^M + '^{ai,a2,a3}- 

This suggests that the order in which modules are taken into account plays 
an important role in the definition of the composition operator. Thus, we 
will focus on an ordered 'partition vr = (X\, . . . , X^) of A, i.e., a partition 
of A provided with a strict total order and represented by a sequence, 
called a modular organisation, preserving (1). Furthermore, we would like 
to be able to "fold" contiguous modules in tt in order to deal with them 
as with a single module^. As a consequence, we require a modular or- 
ganisation to support folding and to be so the composition operator is 
associative according to the order in tt. 

In order to form a modular organisation, the modules and their order in tt 
should satisfy some conditions related to their dynamics. Intuitively, two 
disjoint sets of agents Xi and Xj, (i < j), can be modules in tt, either if 

^ A partition of a set A is a set of nonempty disjoint subsets of A which covers A. 
The folding of modules corresponds to the union of these modules. 
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they do not regulate each other, or if Xi regulates Xj. In both cases, we can 
remark that the equilibria of Xi should embed the asymptotic evolution of 
Xj, which leads to encompass the equilibria of Xj in the equilibria of X^. 
These conditions are expressed by the modularity relation (M-relation) . 

Definition 5. The M-relation^ C 'P(^) x V{A) is defined as follows: 

Xi ^ Xj ^yS'cS: {^x, o ^x,ux,{S'))^x, Q {^x, o ^x,ux,{S')). 

Proposition 2. Let S' be a subset of S. For all Xi, Xj subsets of A, we 
have the following properties: 

1. Xi ^ Xj ^ ^XiUX,{S') C {^x, o QxiUX,{S')); 

2. Xi ^ Xj ^ ^x,ux,{S') = ^x, o ^x,ux,{S'). 

Proof 1. (^) Let S' C S, s e ^x,ux,{S') and s ^ ^x,{nx,ux,{.S')), 
and s' G !fxi(^XjUXj ({•§})). By definition of equilibrium, s' ^x^ux 
Now, we have: 

- Ms" G (s' ^xj : s" e 'l'x,o^x,uXji{s'}), by definition of equilibria; 

- Vs" G (s' ^xj) ■■ s" G ^Xi o ^XiUXjiis'}), by definition of 
As a consequence, s' y^XiUXj which leads to a contradiction. 

(^) Let S'CS,se ^x,ux,{S'), and ^x,ux,{S') C WxA^x,ux,{S')). 
Then, by hypothesis, we have: 

Vs" G {s -X,) : s" G ^x.ux,{S') A s" G ^xA^x,ux,{S')), 
which means that s, s" G <Pxi o ^XiUXjiS'). Thus, Xi Xj. 

2. Prom Proposition 2.1 and since ^x,ux,iS') = ^^x,ux, (l^x.ux, (5"))- □ 

In this context, a modular organisation can be defined as: 

Definition 6. A modular organisation is defined as an ordered partition 
{Xi, . . . , Xm) of A such that, for alll < i < m: (U}=i -^j) ^ -^i- 

Prom Definition 6, Proposition 3 below states that being a modular or- 
ganisation is preserved by any folding of its contiguous parts. 

Proposition 3. Let tt = {Xi, . . . ,X^) be a modular organisation. For 
all 1 < i < j < m, {Xi, . . . , IJ'^^^ Xk, Xj^i, . . . , X^) is a modular 
organisation. 

* Notice that the Af-relation is reflexive but, in general, not transitive and that, when 
its definition holds for some S' C 5*, it holds for any subset of S' . 
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Proof. Let vr = (Xi, . . . , Xj, Xj+i, . . . , Xm) be a modular organi- 

sation. Let X = IJA;=i^fc- We want to show that (Xi, . . . , U 
Xj+i, . . . ^Xsfn) is a modular organisation. By definition 6, we have: 

{X U Xi) ^ (2) 

and: 

X Xi. (3) 

We want to show that: (2) A (3) =^ X {Xi U Xi+i). First, by 
Proposition 2.2, we can write: 

^xuXi o ^xuXiUXi+i = ^xuXiUXi+i by (2), (4) 
^x o ^xuXi = ^xuXi by (3) . (5) 

Thus: 

'^x ° '^xuXiUXi+i = li^x o (!^xux, o 'f'xuXiUXi+i) by (4) 

= !^'xux, o !^xux,ux,+i by (5) 
= i^'xuXiUXi+i by (4), 

which is the expected result. From Proposition 2.2, wc can deduce that: 
X [Xi U Xj+i). Iteratively, we show that X [J-^^,. X^. As a result, 
{Xi, . . . , Xi-i,\Jl,^- Xk, Xj^i, . . . , Xm) is a modular organisation. □ 

In the literature [5,9,10], modules are frequently related assimilated to 
SCCs of interaction networks. Although these works focus on structural 
arguments only, it turns out that they are compatible with Definition 6. 
Indeed, any topological ordering of SCCs is actually a modular organ- 
isation. For instance, ({ai}, {02, as}, {04}) is a modular organisation of 
the interaction network presented in Figure 1. Similarly, any other struc- 
tural or dynamical property could be helpful in the research of modular 
organisations. In what follows, we present an approach addressing for- 
mally this aspect. As a result, we show that, in particular, the structural 
decomposition in SCCs makes sense and may be improved by a deeper 
analysis leading to the decomposition of SCCs in elementary modules (see 
Section 4). 
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3.2 Regulation and modularity relation 

In this section, we focus on the relation between the notions of regulation 
and M-relation. More precisely, and unexpectedly, we will see in Lemma 1 
that this relation refers more directly to the absence of regulation than to 
the regulation itself. 

Lemma 1. For any Xi, Xj subsets of A: -'{Xj — > X^) =^ X^ -w Xj. 

Proof. By Definition 4, for any Xi,Xj C A and for any s,s' G S, we 
have ^{Xj Xi) A (s ^a\Xj s') =^ {s ^xj {s' ^xj- This 
property is obviously preserved at equilibria. Indeed, for any s, s' G S, we 
have -.{Xj Xi) A (s s') =^ "^x^Cs')- Thus, the 

restrictions ifxj(s) and l^Xjls') to Xi are identical. Then, the evolution by 
Xj from the equilibria of Xi remains in the equilibria of Xj. Hence, we 
get that -.{Xj — > Xi) =^ Xi-^ Xj. □ 

Also, as a remarkable fact, according to Definition 6, Theorem 1 provides 
a connection between structural properties of a regulation graph and the 
corresponding modular organisations (possibly reduced to a single mod- 
ule). 

Theorem 1. Any topological ordering of the SCC quotient graph of a 
regulation graph is a modular organisation. 

Proof. Observe that, in the SCC quotient graph G of a regulation graph, 
Xi — > Xj always implies that ^{Xj — > Xi), because of the acyclicity of 
G. Thus, folding contiguous modules with respect to any topological order- 
ing preserves the absence of regulation. As a consequence, if {Xi, . . . , X^) 
is a topological ordering of G, for all i,j G IN such that 1 < i < j < m, 
we have ^{Xj — > Xi), and by Lemma 1, Xi Xj. □ 

3.3 Composition operator 

In this section, we present the successive steps leading to the definition of 
the composition operator 0. From (1), is a binary operator that applies 
on the equilibria of parts Xi and Xj of vr, with i < j. Thus, its definition 
is based on the attractors of ^Xi which correspond to terminal nodes of 
the SCC quotient graph of ^Xi, namely G/^x- ' where is the equiv- 
alence relation on states defined as s ^Xi s' = {s ^^st^ ^ i^' ^Xj 
For any S' Q S, an attractor of ^x^ complies to: [sj^j^. C ^Xi{S') (see 
Figure 4. a). Moreover, for all S' C S, we denote by: 
— [5']^^, = {[s]^-^. \s G 5"} the sets of equivalence classes of ^Xj in S'; 
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(a) (6) 
Fig. 4. Successive steps leading to the definition of the composition operator 0. 



- [sUx, [-^x,U^^ [s%^^ = 3s G G [s%^^ : s ^x, s' the 

evolution by agents of Xj on these equivalence classes. 
We define an operator [•f'x that computes the set of equilibria of 
[-^Xj\^X' ™ Q l^x- (see Figure 4.b) as follows: 

Operator is thus defined as: 

"^Xi = Flatten o o ^x^ , (6) 

where, for any set E C V{S), Flatten(i?) = IJeeE ^- Thus, if applied to a 
set of attractors. Flatten gives the underlying set of equilibria^. As a result, 

it can easily be seen that does compute the set of states belonging to 
the attractors of Xi which are also the equilibria of ^Xj ■ 

Lemma 2 below shows that the global equilibria of X^ U Xj are obtained 
using the composition ^Xi (Zi ^Xj ■ 

Lemma 2. For all Xi, Xj disjoint subsets of A: 

Xi ^ Xj ^ ys' CS: ^x,ux,{S') = i^x, 0^Xj)o f2x,ux,{S'). 

Proof. Let Xi,Xj C A such that Xi ^ Xj and S' a subset of S. 
(C) First, let us show ^x,uXj{S') Q (fx, ^Xj) o ^x,uXj{S'). From 
Proposition 2.2, we know that 'I^Xi o ^XiUXjiS') = ^XiUXjiS'). Thus, 
Vs G ^x^uXjiS'), we have [s]^^^, C 'l^x^uXj{S'). Similarly, an evolution 

" For example, Flatten({{00, 01}, {10}}) = {00,01,10}. 
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by Xi U Xj from an attractor of Xi remain in the same attractor except 
potentially with evolutions by Xj. We have then, for all s, s' G ^xiuXj {S')- 

[sUx,[-^XiUXj]UxyUx, ^ 

Now, since both s and s' belong to attractors of X^ U Xj (by hypothesis), 
if there exists an evolution by Xj from [s]^^. to [s']^^ , there exists 
obviously another path labelled by Xj from [s']^-^. to [s]^-^. . Hence, for 
all s,s' e ^XiUXjiS'), we have: 

{[SU,^ PX,]:.,^ [S'U., ) = {[SU.^ [--X,UX,]U,^ [S'U.^ ) ^ 

As a result, we have [s]-jf. G [^Xj]^xS^^XiUXj{S')]^xJ^ fo^' s G 
^XiUx iS')- Moreover, since from Proposition 1, operator ['Fx ■]^x i^') 
monotonous and since '^XiUXjiS') ^ ^XiUXj{S'), for all s G ^xuriS'), 
we can write that [s]^^, G [^Xjj^xX^XiUXjiS')). Now, since s G [sj^xo 

Vs G ^XiUXjiS'), we have s G Flatten o [^]^^. (^^x^ux,- (5'))- Prom (6) 
and Proposition 2.2, we can write: 

^Xi ifx,- o f2xiUX^{S') = Flatten o [^]^^^ o ip-x, o i7x,ux,(.S') 

= Flatten o [^,Ux, ° ^x^ux.iS'). 

Hence, for all s G 'l^xuriS'), wc have: s G l^x '^y ° ■^^xuyl-S")) which 
corresponds to the following inclusion: 

^xuy{S') C^x0^yo nxuY{S'). 



(D) Now, let us show {'Fx, Fxj) o ^x^uXjiS') C 'Fx.uXjiS'). To do 
so, let us consider a state s G (Fx^ ^Xj) o ^^x^uXj (5")- From (6) and 
Proposition 2.2, we have [s]^^. G [!fxj^^. o [Vx^ o nXiUXj{S')U^.. 
Now, consider [[s]^^ . By definition of attractors, for all si,S2 G 
Flatten o Flatten we have si ^XiUXj This means that 
= ux '^rid, as a consequence, that s G '^xuviS')- As a 

result, the inclusion (Fx Fy) o fixuviS') Q ^xuy{S') holds. □ 

As a main result. Theorem 2 shows how can be used on a modular 
organisation vr to obtain a modular computation of global equilibria of an 
interaction network. 
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Theorem 2. Let A' = UI=i Xi A he a set of agents such that {Xi, . . . , 
Xm) is a modular organisation, we have: 

^A' = i^x, ... ^xj o Qj,,. 

Proof. This proof is made directly, using Definition 6 and Lemma 2. Since 
TT = {X\, . . . , Xm) is a modular organisation, it is folding preserving and 
Ui^l^ Xi Xjn. Then, using (6) and Lemma 2, we have: 

= Flatten o [<Fx^]^ , o tZ^, ,™ i ^ o Oa' 

Ui — l '2 = 1 » 

= Flatten o o ^^{JTJ^" ® '^^—i) ° -^UfcT'^i ° 

= Flatten o o <Z^x^_ J o ^2^, 

= Flatten o [^U o [f^.U o ^^^-2 o 1?^, 

Ui=i Ui=i '-'1=1 » 

= Flatten o [^U 0...0 [^]^ o ^x, o Qa'- 

Ui=i 

As a result, we obtain: 

^A'{S') = {^x,0..-0^xjo QA'iS'), 
which is the expected result. □ 

4 Elementary modular organisation 

Informally, a module is elementary if it is not separable, i.e., if the equilib- 
ria of each of its agents depend entirely on the equilibria of all the others. 
For instance, consider negative circuits that lead to asymptotic sustained 
oscillations [19]. In such regulation patterns, the equilibria of an agent 
cannot be encompassed into that of the others because, in order to reach 
its own equilibria, each agent evolves from the equilibria of all the others. 

In this context, a modular organisation provided by some topological or- 
dering of the sec quotient graph (see Theorem 1) does not always provide 
an elementary decomposition. So, it may be not satisfactory since modules 
often support specific biological functions, revealing elementary modules 
appears important, not to say essential, to highlight basic biological func- 
tional mechanisms. 
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Fig. 5. A separable interaction network with modular organisation ({6},{o, c}). 



Figure 5 depicts an interaction network rj composed of three agents ai, 
02 and as with the associated strongly connected regulation graph. Its 
underlying state graph shows that the global dynamics of rj leads to two 
attractors, stable state {111} and limit set {000, 100, 101,001}. It is easy 
to see that {02} {^17^3} and that this Af-rclation is (obviously) pre- 
served by folding, because there are only two modules. Hence, ordered 
partition ({02}, {oi, 03}) is a modular organisation of 77. 

However, separability is not possible in general as illustrated in Figure 6. 
Indeed, starting from modular organisation vr = ({0.1}, {027 ^3}) obtained 
from the SCCs, the separation of {02, 03} should lead to the ordered parti- 
tion tt' = ({ai}, {02}, {o^a})- The condition for w' to be a modular organi- 
sation is that {fli, 02} {03}) ^-6-) the evolution by 03 from the equilibria 
of {oi, 02} have to be included in the equilibria of {oi, 02}. We can observe 
that the attractors for {01,02} are {010,011} and {100,101}, while the 
evolution by 03 from either Oil or 101 leaves the attractors of {01,02}, 
which means that vr' is not a modular organisation. As a consequence, 
{02, 03} cannot be separated. 

Hence, the separation condition of a module in tt is not local to this 
module but depends on the module "context", that is the global equilibria 
{i.e., '^yi-i modules that precede it in tt. Deciding the separabil- 

ity of Xi into Xl and Xf implies checking two conditions: IJfc=\ -^k ^1 
andUfe=i^fcUXi ^Xf.Of course, the complexity of the underlying com- 
putation is exponential in the size of tt and depends also on the position 
of Xi in TT. Nevertheless, brute-force computation may be used in prac- 
tice for small interaction networks (of about 15 agents). A more efficient 
method allowing to go beyond this limitation is, for the moment, an open 
question. 
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no.2{s) = (Sai A Sas) V (Sai A "'Saa) V (Saj A -^Sa^) 
Vaais) = (Soi A Saa) V (-iSaj A Sag) V (Saa A Sag) 




Fig. 6. Example of non-separability of {02,03} in tt = ({ai}, {02,03}). 

5 Conclusion 

We developed a formal framework for the analysis of the modularity in 
interaction networks assuming asymptotic dynamics of modules and en- 
abling their composition. We exhibited modularity conditions governing 
the composition of modules and an efficient computation method such that 
the global equilibria of interaction networks are obtained from the local 
ones. Moreover, we confirmed that usual assumptions identifying modules 
to SCCs have a strong motivation coming from theory. Of course, the 
next step should be finding ingenious algorithms for searching elementary 
modular organisations. Then, since this work provides a rigorous setting 
for studying other questions around modularity, it would be of interest 
to exhibit elementary biological components, which would be relevant, in 
particular, for synthetic biology [22]. Also, questions related to robust- 
ness and evolution could be tackled thanks to the modular knowledge of 
interaction networks. 
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